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Motivated by the gravity/fluid correspondence, we introduce a new method for characterizing nonlinear grav¬ 
itational interactions. Namely we map the nonlinear perturbative form of the Einstein equation to the equations 
of motion of a collection of nonlinearly-coupled harmonic oscillators. These oscillators correspond to the quasi¬ 
normal or normal modes of the background spacetime. We demonstrate the mechanics and the utility of this 
formalism within the context of perturbed asymptotically anti-de Sitter black brane spacetimes. We confirm in 
this case that the boundary fluid dynamics are equivalent to those of the hydrodynamic quasinormal modes of 
the bulk spacetime. We expect this formalism to remain valid in more general spacetimes, including those with¬ 
out a fluid dual. In other words, although borne out of the gravity/fluid correspondence, the formalism is fully 
independent and it has a much wider range of applicability. In particular, as this formalism inspires an espe¬ 
cially transparent physical intuition, we expect its introduction to simplify the often highly technical analytical 
exploration of nonlinear gravitational dynamics. 

PACS numbers: 04.70.Bw 

I. INTRODUCTION 

Can spacetimes become turbulent? Direct numerical sim¬ 
ulations of large asymptotically anti-de Sitter (AdS) black 
holes [1] and their holographically dual fluids [2, 3] have pro¬ 
vided convincing evidence that this is the case. This phe¬ 
nomenon, perhaps counterintuitive at first glance,* can be un¬ 
derstood through the gravity/fluid correspondence [4-6]. This 
correspondence links the behavior of long-wavelength pertur¬ 
bations of black holes in AdS to viscous relativistic hydro¬ 
dynamics, and its regime of applicability can include cases 
of high Reynolds number on the fluid side. Spacetime turbu¬ 
lence then follows from turbulence in the dual fluid [2, 6]. On 
the gravity side, a high Reynolds number corresponds to dis¬ 
sipation of gravitational perturbations that is weak when com¬ 
pared with nonlinear interactions. It is therefore not surprising 
that it arises in the vicinity of asymptotically AdS black holes, 
which can have relatively long lived quasinormal modes. 

The observation of gravitational turbulence in AdS moti¬ 
vates a further question: Can one analyze this striking non¬ 
linear behavior directly in general relativity without relying 
on the existence of a holographic dual! That is, rather than 
borrowing from the dual hydrodynamic description—and any 
restricted regime of applicability—can one establish a bona- 
fide description of turbulence as a perturbative solution of 
the Einstein equation? Recall that turbulence is a nonlinear 
phenomenon characterized, in particular, by cascades of en¬ 
ergy (and sometimes enstrophy) between wave numbers. It is 
therefore delicate to fully capture this behavior within ordi¬ 
nary perturbation theory without carrying it out to sufficiently 
high orders and performing a suitable resummation [3]. In 


' Due to a crucial difference: the Einstein equation is linearly degenerate as 
opposed to truly nonlinear as is the case of e.g., the Navier-Stokes equa¬ 
tions. 


order to take into account the essential gravitational self¬ 
interactions of perturbations that are present in the Einstein 
equation we will require a more general perturbative frame¬ 
work. 

In this work we introduce a nonlinear coupled-oscillator 
model to describe the interaction of quasinormal or normal 
modes of a background spacetime, in particular their mode¬ 
mode couplings. This proposal is a natural generalization of 
our earlier study of nonlinear scalar wave generation around 
rapidly-spinning asymptotically flat black holes [7], where the 
back-reaction on the driving mode was neglected (we account 
for it properly in this paper). This previous model illustrated 
that the onset of turbulence in gravity does not require the 
spacetime to be asymptotically anti-de Sitter^. In the non¬ 
linear oscillator model presented here, the coupling between 
modes is accounted for explicitly and in real time as opposed 
to implicitly through a recursive scheme. Therefore the equa¬ 
tions of motion provide solutions that are valid over longer 
time scales. 

Within this model, nonlinear gravitational perturbations are 
described by excitations of modes (quasinormal or normal). 
Eor a given background spacetime, the collection of modes 
is parametrized by a particular set of frequencies, damping 
rates, and, at the nonlinear level, mode coupling coefficients. 
Through these parameters, we can quantitatively compare and 
contrast signatures of nonlinear gravitational perturbations in 
different backgrounds, in the same way that frequencies and 
damping rates alone characterize linear perturbations. In this 
way we can gain a better understanding of nonlinear interac¬ 
tions and associated phenomena (such as turbulence) in gen¬ 
eral relativity. The route taken when constructing this for¬ 
malism essentially offers a new perspective on how to deal 


^ In analogy to hydrodynamics, it is of course necessary to be in the regime 
of high gravitational Reynolds number. 
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with nonlinear metric perturbations that is conducive to intu¬ 
ition building. This compares favorably with more traditional 
methods, where one has to contend with difficult technical de¬ 
tails that often mask the underlying physics. 

To provide a concrete example, we will apply our meth¬ 
ods to study nonlinear perturbations of an asymptotically AdS 
black brane. The gravity/fluid correspondence applies in this 
case and the resulting coupled-oscillator system may be com¬ 
pared against the dual fluid. We hnd that our equations are 
consistent with the relativistic hydrodynamic equations pro¬ 
vided by the duality. Although the agreement is expected, our 
derivation provides an explicit demonstration and a natural 
physical interpretation of the observed phenomena in terms 
of quasinormal modes. We emphasize that the derivations in 
the gravity and fluid sides are independent of each other, and 
so the treatment for gravitational perturbations does not de¬ 
pend on the existence of a dual fluid and can be applied to 
more general spacetimes. 

In the interest of caution, we recall that quasinormal modes 
do not form a basis for generic metric perturbations (see [8] 
for a recent discussion). For instance, consider linear pertur¬ 
bations of the (asymptotically flat) Kerr spacetime as an ex¬ 
ample (see also discussions in Sec. II). The signal sourced by 
some matter distribution comprises quasinormal modes, the 
late-time “tail” term, as well as a prompt piece that travels 
along the light cone. In this sense, our formalism is approx¬ 
imate as we consider only the quasinormal mode contribu¬ 
tions. However, in many cases, such as the ringdown stage 
of binary black hole mergers or when considering long wave¬ 
length perturbations of an asymptotically AdS black brane, it 
is sufficient to track only the quasinormal modes, as they are 
the dominant part of the signal (see, e.g., [9], for a related dis¬ 
cussion). In more general scenarios, we can always check the 
validity of our approximation by estimating the magnitudes of 
the other contributions. 

This paper is organized as follows. In Sec. II, we intro¬ 
duce the general formalism of the nonlinear coupled-oscillator 
model, and compare it with traditional methods for handling 
nonlinear gravitational perturbations. In Sec. Ill, we briefly 
review the asymptotically AdS black brane spacetime and the 
gravity/fluid correspondence, and we analyze the boundary 
fluid in the mode-expansion picture. In Sec. IV, we apply the 
general formalism to the specific case of the asymptotically 
AdS black brane. We conclude in Sec. V. The gravitational 
constant G and the speed of light c are both set to one, unless 
otherwise specified. Appendices are provided to elaborate on 
certain details. 


II. GENERAL EORMALISM 


In this section, we begin by reviewing the traditional ap¬ 
proach to solving the Einstein equation using ordinary per¬ 
turbation theory and assuming a series expansion in the per¬ 
turbation amplitude. This method might not lend itself to 
easily capturing relevant phenomena like turbulence. In the 
case where the linearized dynamics take the form of indepen¬ 
dently evolving normal or quasinormal modes (in the absence 


or presence of dissipation, respectively), we then show how 
the nonlinear Einstein equation can be represented as a set 
of coupled oscillator equations, which is analogous to treat¬ 
ments of the Navier-Stokes equation in fluid dynamics, and is 
indeed capable of cleanly capturing turbulence. Eor simplic¬ 
ity, we restrict our discussion to vacuum spacetimes, but it is 
straightforward to generalize the analysis to spacetimes with 
a cosmological constant. 


A. Ordinary perturbation theory 

Given any metric one can split it into the sum of a 
“background” metric and a “perturbation”, 

g^v = g^v ■ ( 2 . 1 ) 

Without invoking any approximation, the vacuum Einstein 
equation may then be written as 

oo 

(/ » + R^i^l (g^ , h) -f ^ W (gi", h) = 0, 

71=3 

( 2 . 2 ) 

where R^^v{g^,h) denotes the nth order Ricci tensor expanded 
about g®y. Explicitly, the linearized and second order terms 
are 
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In these expressions, covariant derivatives associated to the 
background metric g®y are denoted by vertical lines. The 
background metric is also used to raise and lower indices. 

As described in [10], ordinary perturbation theory assumes 
the existence of a one-parameter family of solutions g/^y(e), 
where gfiv(O) = g^y, and /r^y depends differentiably on e. 
One can then Taylor expand the perturbation, 

h^y = ehl!l + e^h^^l + --- . (2.5) 

Perturbative equations of motion of order n follow by differ¬ 
entiating the Einstein equation (2.2) n times with respect to e, 
and then setting £ = 0. At zeroth order we have simply 

( 2 . 6 ) 


so that is a vacuum solution itself. 

At first order in £ we have the linearized Einstein equation. 


(2.7) 
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It is generally much easier to solve this equation (after making 
appropriate gauge choices and imposing boundary and initial 
conditions) than it is to solve the full Einstein equation. Then 
for sufficiently small e, + e/jj/v should be a good approx¬ 
imation to g^v{e). 

This procedure may be continued to higher orders. For in¬ 
stance, at second order, we obtain 

R^,!kg^.h^^^) = -Rfkg^,hW). ( 2 . 8 ) 

The second order perturbation is seen to evolve in the back¬ 
ground spacetime and it is sourced by the first order so¬ 
lution 

Generically, this approach reduces the nonlinear problem to 
a series of linear inhomogeneous problems of the form 

4v) = Sgl (g®;/r('),...,• (2.9) 

Thus, at each order, one solves a linear partial differential 
equation with a source, subject to appropriate boundary con¬ 
ditions and gauge choices. The left hand side of the equa¬ 
tion at order n consists always of the nth order perturba¬ 
tion evolving linearly in the background spacetime g^y. 
The source term involves only already-solved lower or¬ 
der pieces h^v for < n, so a higher order perturbation does 
not backreact on one of lower order. Moreover, since the 
nth order perturbation evolves in the zeroth order background 
metric—not the (n — l)th order metric—the efficient capture 
of parametric resonance type effects is precluded [3, 7]. (Of 
course, with enough intuition, it may be possible to identify 
this behavior through a suitable resummation of perturbations 
of sufficiently high order.) In following this program, the cal¬ 
culations are quite involved and the gauge choices at different 
orders are often subtle (see e.g., [11-14]). In the specific con¬ 
text of extreme mass ratio binaries, recent examples of this 
program are given in [15, 16]. 


B. Larger perturbations 

After iterating the above procedure to any given order, the 
resulting perturbative metric should be a good approximation 
to gab{£) for sufficiently small £. However, in certain situa¬ 
tions one may be interested in studying systems with larger 
(but still small) values of e, where the Taylor expansion (2.5) 
either fails to converge or would require a large number of 
terms to obtain a good solution. Typically the perturbative so¬ 
lution would be valid for a short time, but for long times secu¬ 
lar terms might dominate. Therefore, a more suitable scheme 
would be required. In, for example, the context of the Navier- 
Stokes equation, ordinary perturbation theory might be capa¬ 
ble of capturing the initial onset of turbulence, but it would be 
ineffective in capturing fully developed turbulence (and like¬ 
wise for gravitational turbulence [3, 7]). 

In order to characterize the nonlinear dynamics in gen¬ 
eral relativity in a more efficient and transparent manner, we 


present here an alternative way of obtaining approximate solu¬ 
tions that is better suited for exploring certain nonlinear phe¬ 
nomena such as wave interactions and turbulence. We assume 
as before that satisfies the vacuum Einstein equation. But 
then, rather than Taylor expanding /i^v as in (2.5), we consider 
the full metric perturbation h^y, and we attempt to solve di¬ 
rectly a truncated version of (2.2). In fact truncation at second 
order, 

( 2 . 10 ) 

captures the essential nonlinearities of interest to us here. We 
note that our formalism could straightforwardly be extended 
to higher orders, but for simplicity we restrict to second order 
nonlinearities here. 

To summarize, instead of solving a tower of inhomoge¬ 
neous linear equations (2.9) we solve a nonlinear equation, 
but we neglect the higher order nonlinearities. Instead of deal¬ 
ing with gauge issues at each order, we have only to impose 
the gauge condition once on hf^y. Of course, the truncation of 
the Ricci tensor is not a tensor itself so the equation (2.10) is 
not gauge invariant. But it should be sufficient to the order we 
are working {0{h^)). As we shall see, this approach readily 
captures the nonlinear mode coupling effects of interest to us. 

In general it will be very difficult to solve (2.10), even 
neglecting the higher order nonlinearities as we have done. 
However, as we describe in the following subsection, in cases 
where the linear dynamics is dominated by the evolution of 
normal or quasinormal modes, (2.10) reduces to a system of 
nonlinearly coupled (and possibly damped) oscillators. 

C. Expansion into modes 

We now restrict consideration to background spacetimes 
whose linear perturbations are characterized (for some region 
of spacetime) by a set of modes (normal or quasinormal). In 
this case the first order metric perturbation may be written 

hk{t,x) ^ Y.kj{t)3fk^\x)+qfit)3^lii^\x)], ( 2 . 11 ) 

j 

with 

qj it) = Aje-‘^ffi qjit) = Bje'^’j'. ( 2 . 12 ) 

The background spacetime is assumed to be stationary and the 
t coordinate is the associated Killing parameter. Modes al¬ 
ways occur in pairs with frequencies (Oj and —cOj, so we have 
organized the summation above along these lines, labeling 
each pair with a multi-index j (denoting both the transverse 
harmonic and radial overtone). The associated spatial wave 
functions are denoted 3fkv^\x). Finally, q^ and {Aj,Bj} are 
the displacements and the amplitudes for modes y±, respec¬ 
tively. As hf^y must be real at all time, we expect that {Aj,Bj} 
(as well as are conjugate to each other. 

The reason we organize our modes into pairs in (2.11) is 
to emphasize that all modes must be included in the nonlin¬ 
ear analysis; many linear analyses use symmetry arguments to 
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only treat modes with 5R(a)) > 0 [17]. In the case of normal 
modes, the mode functions are degenerate and (Oj G K, 

so we take qj = qj For quasinor¬ 
mal modes, the radial dependence of along with the 

dissipative boundary conditions at the horizon and/or infinity, 
fixes the time dependence of the mode uniquely. Any “degen¬ 
erate” mode in this case must therefore have (Oj = — so 
the frequency is purely imaginary, and the multi-index j de¬ 
scribes just a single mode. We analyze these cases separately 
from the non-degenerate case in the following sections. 

Frequencies of quasinormal modes have nonzero positive 
imaginary part, which implies an exponential time decay as 
a result of energy dissipation. In addition, this complex fre¬ 
quency means that the mode functions generally blow up at 
spatial infinity and the horizon bifurcation surface. How¬ 
ever, as physical observers effectively lie near null infinity, 
the quasinormal-mode signals they observe are finite and the 
modes are indeed physical perturbations of the spacetime. For 
such observers, the sum in (2.11) can become a good ap¬ 
proximation over finite time intervals, although we remind 
the reader that quasinormal modes do not form a complete 
basis for generic metric perturbations^. Additional contribu¬ 
tions to the metric can arise at late times from waves being 
scattered by the background potential at large distances (the 
“tail” term), or at early times from a prompt signal (on the 
light cone) from the source (see, e.g., [17-20]); we collect 
these into the “residual part”. 

In this paper our focus is on mode-mode interactions and 
the associated coupling coefficients. We will therefore not 
consider the nonlinear interactions between the modes and 
the tail and prompt components of the metric perturbation. 
We caution, however, that such couplings need not always be 
small. While they are small for perturbations of AdS black 
branes in the hydrodynamic limit (which we analyze below), 
readers should keep in mind that they will lead to additional 
contributions to, e.g., Eq. (2.15) below. Furthermore, ques¬ 
tions as to how quasinormal modes are excited by moving 
matter, or how to compute the excitation factors for these 
modes based on some arbitrary initial data are also beyond 
the scope of this work (see [19, 21, 22] and Appendix D). 

With these observations in mind, following the discussion 
in Sec. IIB we write the full metric perturbation as 

Vv(l,x) = Y^[qj{t)3f^t\x)+qj{t)^^i^\x)] 
i 

-f residual part, (2.13) 

but now generalizing the coefficients Aj and Bj to be functions 
of time, 

qf{t) =Aft)e-'<^A^ q+{t) = Bj[t)e‘^V . (2.14) 

Our task is to determine the nonlinear evolution of quasi¬ 
normal modes; in other words, to evaluate the time depen¬ 
dence of qJ. Addressing this task is generally nontrivial as 


^ This qualification is represented by the use of the notation in (2.11) 
(see, e.g., [18]). 


it requires the proper separation of the quasinormal modes 
from the residual part of the full metric perturbation. For 
Schwarzschild and Kerr spacetimes this is achievable by in¬ 
voking the Green’s function technique (Appendix D), whereas 
the generalization of this approach to generic spacetimes re¬ 
mains an open problem. To present the coupled-oscillator 
model, we apply an alternative strategy of plugging (2.13) into 
the truncated Einstein equation (2.10) and projecting our the 
spatial dependencies, thereby obtaining mode evolution equa¬ 
tions. This method is most accurate for dealing with normal¬ 
mode evolutions and cases where the residual parts are neg¬ 
ligible (for example, see Sec. IV). In more general scenarios, 
we shall make several additional approximations (such as ne¬ 
glecting certain time derivatives, neglecting the residual part) 
to single out the ordinary differential equations for q^. We 
also caution that since the set of modes generally does not 
form a complete basis, the resulting h^v is still only an ap¬ 
proximate solution to the truncated Einstein equation. Eor 
simplicity, hereafter we shall not explicitly write down the 
residual part in the equations. 

Upon substitution, the truncated Einstein equation (2.10) 
takes the form 

E E [pj (*)^/] 

j ■'-i 



Here pj, rj, and cjj are tensor functions of the spatial coor¬ 
dinates, and they depend on the background metric as well as 
the corresponding wave function of the quasinormal mode. 
The right hand side of the equation has a complicated x- 
dependence that we have suppressed. 

We would now like to project Eq. (2.15) onto individual 
modes to obtain equations for a set of nonlinearly coupled 
oscillators in the form of 




= S‘ 


4qf) > ( 2 - 16 ) 


for each j and s. In order to do so we require a suitable set of 
projectors. If, along any of the dimensions transverse to the 
radial direction, the background metric possesses a suitable 
isometry group so that this part of the wave function is de¬ 
scribed by tensor harmonics (Eourier modes, tensor spherical 
harmonics, etc.) then it is easy to project out this part by using 
an inner product. The remaining part (generally including the 
radial direction) is however more problematic. 

It is often the case that the equations can be written in the 
form of a standard eigenvalue problem, ‘F = —A'P. Eor nor¬ 
mal modes, one can define an inner product {xIb) with re¬ 
spect to which A is self-adjoint, and the modes are orthogonal. 
One can then use this inner product to define the projector. Eor 
dissipative systems with quasinormal modes, the eigenvalues 
are complex and A cannot be self-adjoint. Another problem 
is that often | \ —oo at the dissipative boundaries of the 

system. Nevertheless, it is still possible to define a suitable 
bilinear form, with respect to which A is symmetric [7, 23- 
28]. This bilinear form involves an integral of XB without any 
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complex conjugation so symmetry of A does not imply that the 
eigenvalues are real. Furthermore it is still necessary to appro¬ 
priately regulate the integration to eliminate divergences. The 
bilinear form may be regarded as a “generalized” inner prod¬ 
uct, and be used as such. In particular, it may then be shown 
that (J"-'* I = 0 for cOj ^ (Ok, and this orthogonality leads 
to a suitable projector. 

In the general case (such as the coordinate system we use 
in Sec. IV) it is not necessarily possible to re-write the equa¬ 
tion as a standard eigenvalue problem. Nevertheless, we can 
still define a generalized inner product and use it to project 
the equation onto modes. It may be that the modes are not 
orthogonal with respect to this inner product, in which case 
the projection of the left hand side of ( 2 . 10 ) contains contri¬ 
butions from additional modes beyond the desired projection 
mode. After performing projections onto all modes, it would 
then be necessary to diagonalize the system to obtain a set 
of equations of the form (2.16). This is possible by apply¬ 
ing procedures described in Sec. IIC 1 to remove “unphysical 
modes” and reduce the order of the differential equations. At 
this point, it is worth noting that in principle any inner product 
which leaves this set of equations non-degenerate fits our pur¬ 
pose. However, in order to minimize the error from neglecting 
the residual part, it is good practice to adopt an inner-product 
suitable for eigenvalue perturbation analysis (see Sec. IV B for 
a concrete example of such an inner-product). 

With the equations decoupled as in (2.16) with a suitable 
generalized inner product, we can now substitute in Eq. (2.14) 
for qj. We obtain, 

a-^Aj + b-jAj = SJ iAk,B ,), (2.17) 

a+Bj+bpj = S+iAk,B,), (2.18) 

where = ^7 —'2.1(0 and St = +2i(0*a^. We have 

used the fact that and e'“j‘ are homogeneous solutions 

to simplify the left hand sides. The “source” terms on the 
right hand sides are quadratic in A^ and Bk- We have dropped 
quadratic terms involving derivatives of A^ and Bk in as we 
expect them to be smaller than quadratic terms not involving 
derivatives. Indeed Eqs. (2.17)-(2.18) already indicate that 
time derivatives of the coefficients are of quadratic order in the 
perturbation amplitudes, so that, e.g., terms on the right hand 
side of the form A<;A/ would be of cubic order. In general, the 
nonlinear terms will then be of the form 


1. Non-degenerate modes 

The non-degenerate case applies to quasinormal modes 
only. We immediately see from examining (2.17)-(2.18) that 
with = 0, {Aj,Bj} = constants are solutions. This is by de¬ 
sign as (2.12) are solutions to the linearized equations. How¬ 
ever, if aj 7 ^ 0 the left hand sides of (2.17)-(2.18) are second 
order in time, so that there are additional homogeneous solu¬ 
tions, 

A — ^7 1 1Cl\ T, — ^ i t j n : /r\ .Af\\ 

Ajoc e Bj°^e J J , ( 2 . 20 ) 

which give rise to 

qj oc ^ J t , qj i ■>'}'. (2.21) 

These solutions are clearly not quasinormal modes since when 
combined with the spatial wavefunctions, they do not satisfy 
the appropriate dissipative boundary conditions. In addition, 
if we multiply them with the wave function the origi¬ 
nal linearized Einstein equation is not necessarily satisfied (if 
fly 7 ^ 0 and b^j 7 ^ 0). At the linear level, one can require Aj,Bj 
to be constants to remove these spurious modes. At the non¬ 
linear level, we need a systematic strategy to eliminate this 
extra unphysical degree of freedom. 

Let us first assume that aj 7 ^ 0. Eor clarity we only con¬ 
sider the s = -f modes, but the analysis carries over directly 
to i = —. We will argue that the second time derivative terms 
in equations (2.17) and (2.18) should be dropped. To arrive 
at an intuition for this, first note that we are considering the 
problem of mode excitation in the presence of sources. In 
equations (2.17) and (2.18), the source terms come from non¬ 
linear couplings, but it is more instructive to move beyond 
this particular specialization and consider generic sources. If 
a delta-function source S = ~Xq) is introduced to the 

spacetime, it gives rise to a finite-value discontinuity of the 
quasinormal mode amplitude at f = fo, after which quasinor¬ 
mal modes evolve freely and Aj remains constant (see the ex¬ 
ample in Appendix D). In other words, only Aj jumps at the 
delta source while A,- is unaffected (otherwise it will not re¬ 
main constant in the ensuing free-evolution), so that only Aj 
is needed in a sourced mode evolution equation to account for 
the influence of that source, while A does not in fact contribute 
to the evolution of the physical modes. Eurthermore, dropping 
A also frees us of the unphysical spurious modes, as the evolu¬ 
tion equation is now first order in time. We have subsequently 


Ik ^ 

-hKjPBkBP“k+“i+“j'>' , (2.19) 


— (n) 

where the coefficients KjjJ ' are constants (and similarly for 



We now proceed to separately analyze non-degenerate and 
degenerate modes. 



Mathematically, this physical intuition is reflected in the 
fact that when we integrate (2.17) from fo_ to fo+ with a delta- 
function source at f = fo, we realize that the integration of the 
A term in fact vanishes because A^ (fo_) andA j(fo+) must both 
be zero in order to satisfy the free evolution condition when 
the source vanishes. We note of course that the solutions of 
equation ( 2 . 22 ) no longer strictly satisfy the original equations 
(2.17) or (2.18). However, since both set of equations should 
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be satisfied on physical grounds, A and B terms should be bal¬ 
anced by the residual part of the metric perturbations, which 
is implicit in the left hand sides of (2.17) and (2.18). 

The situation with a* = 0 does not present any of the above 
difficulties as the oscillator equation (2.17) or (2.18) is already 
first order in time, so that 





or 



(2.23) 


In fact, this is the case we shall encounter in Sec. IV when 
we perturb about the anti-de Sitter black brane background in 
ingoing Eddington-Finkelstein coordinates. In that case per¬ 
turbations are described by a first order in time and second 
order in space partial differential equation. 


2 . Degenerate modes 

For a degenerate mode, the two equations in (2.16) for 
s = ± degenerate to a single equation for a, = a7 + Thus 
the 4 degrees of freedom present for a given j that we saw in 
the non-degenerate case reduce to 2 degrees of freedom (or 1 
if Qj = 0). In other words, we do not have any unphysical spu¬ 
rious solutions in the degenerate case, but instead two sets of 
physical solutions with the same spatial wavefunction, which 
should both be kept. The consequence of this observation is 
that in the end, the evolution equation for each mode is of 
first order, and we need not apply the treatment for the A term 
employed in the non-degenerate case. 

Consider first the case where aj ^ 0. As noted earlier, this 
corresponds to a non-dissipative (i.e. normal) mode. An ex¬ 
ample where this occurs is in perturbations about pure anti- 
de Sitter spacetime (without any black hole). (The case of 
coupled scalar field-general relativity perturbations about AdS 
was analyzed as coupled oscillators within the context of a two 
timescale expansion in [29].) 

As discussed before, even for this aj 7 ^ 0 case, the (2.17)- 
(2.18) should reduce to first order, and we show below how 
this is to be achieved. First note that we have 

qj =Aj{t)e-^‘^i^+B,{t)e^^j ‘, (2.24) 

and when we introduced time dependence into Aj and Bj, 
these parameters can in themselves contain and 

factors, so their choices in equation (2.24) are not unique, and 
we have in effect a freedom that we have to fix. The most 
obvious optimal choice is to enforce 

qj = -iO}jAje-'^i’ + ico*Bje^“P. (2.25) 

as a gauge fixing, or equivalently 

Aje-“^J‘ + Bje‘“j‘ = 0 , (2.26) 

which incidentally looks as if we were solving an inhomo¬ 
geneous equation through a variation of parameters method. 
The physical intuition behind this constraint is that Aj and Bj 
change only slowly with time so it is appropriate to regard 


them as “instantaneous” amplitudes. (However, this does not 
constitute a restriction on the solution.) We then have 


Aj = 


(0*qj + iqj 


co + co* 


D „-icop 0) qj - iqj 
ts j = e J - 

^ CO + CO* 


(2.27) 


So far we have not imposed any equation of motion, and after 
substituting in equation (2.16) and walking through the same 
procedure as that presented in Appendix A, we obtain 


^; = 




COj) 




le J ^ 
_ _ c. 

CO*) 


7(® 


(2.28) 


We have thus re-expressed the second order equation (2.16) 
for qj in terms of first order equations for the amplitudes Aj 
and Bj . 

In the case where aj = 0, we have cOj = —wj = —icj/bj, so 
COj is purely imaginary and there is a single degree of freedom. 
There is then no need to distinguish Aj and Bj, so we can set 
Bj = 0. Equation (2.17) easily reduces to 


Sj SjcA^A 


bj bj 


(2.29) 


Equations (2.22), (2.23), (2.28) and (2.29) are our desired 
first order equations of motion. They describe a collection 
of nonlinearly coupled harmonic oscillators. For any suitable 
background spacetime, perturbations are characterized by the 
mode spectrum, the mode-mode coupling coefficients and the 
mode excitation factors. 

Despite being a simplified model in the small amplitude 
limit, the formalism we introduced in this section effectively 
serves as a general platform to quantitatively compare and 
study the nature of nonlinear gravitational phenomena in dif¬ 
ferent spacetimes. A most attractive feature is that the vast 
literature on nonlinear coupled oscillators that has been de¬ 
veloped in other branches of physics can now be applied di¬ 
rectly to the study of gravitational interactions. For example, 
a precursor to the present procedure led to the discovery of the 
parametric instability in the wave generation process in near¬ 
extremal Kerr spacetimes in Ref. [7], which exhibited similar 
properties to the parametric instability in nonlinear driven os¬ 
cillators. In general relativity, another example is furnished by 
the study of perturbed anti-de Sitter spacetimes through a two 
timescale analysis [29] and its connection to the Fermi-Pasta- 
Ulam problem [30, 31]. 

In Sec. IV below (with some details relegated to Appendix 
B), we provide a concrete example on how to implement the 
abstract procedure laid out in this section, using the asymp¬ 
totically AdS spacetime containing a black brane as the back¬ 
ground. The study of this particular case also results in a num¬ 
ber of interesting physical observations, and so has its own 
intrinsic value. For example, we shall see that relativistic hy¬ 
drodynamics admits a similar description to the gravitational 
equations of motion, thus expanding the gravity/fiuid corre¬ 
spondence. Additionally, by connecting it to the fiuid side one 
concludes that the symmetry of K is closely connected to the 







cascading/inverse-cascading behavior in the turbulent regime. 
Hence, this duality mapping provides further evidence and in¬ 
sights for the behavior of turbulence in gravity. 
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III. AdS BLACK BRANE SPACETIMES AND THE 
GRAVITY/FLUID CORRESPONDENCE 


In advance of our analysis of coupled AdS black brane 
quasinormal modes in Sec. IV, here we review the grav¬ 
ity/fluid correspondence and study the black blane perturba¬ 
tions from the fluid side. We first present the background uni¬ 
form AdS black brane solution. We then review the derivative 
expansion method that leads to boundary fluid equations that 
describe long wavelength perturbations. Finally, by Fourier 
transforming the boundary coordinates we re-write the system 
as a set of coupled oscillators to facilitate comparison with our 
later gravitational analysis. For a more complete introduction 
to the gravity/fluid correspondence, interested readers should 
consult the original references [4-6, 32]. 


A. Background metric 


The metric for the d I dimensional uniformly boosted 
AdS black brane is given in ingoing Eddington-Finkelstein 
coordinates by 



—lu^dx^dr + P' 



1 

{brY 


iflUy 


dx^dx^. (3.1) 


where (with = — 1) is some arbitrary constant four 
velocity, r is the radial coordinate and x^ are the boundary 
coordinates. The Hawking temperature of the black brane is 
the constant T = d/{Anb). This metric satisfies the Einstein 
equation 


G;xv+Ag^v=0, (3.2) 

with cosmological constant A = —d{d— l)/2. 

Different choices of correspond simply to different 
Lorentz-boosted boundary frames. In particular, in the case 
where the spatial velocity vanishes, the above metric simpli¬ 
fies to 


d-l 

t/ipj = 2dvdr — Y f {r)dv^ + '"^12 {dx,iY , (3.3) 

i=i 


where /(r) =1 — 1 /{brY ^^d v is the ingoing Eddington- 
Einkelstein coordinate. The horizon is then located at r = \/b. 

If we define the tortoise coordinate r* as dr^ = dr/{Yf{r)) 
and dv = dt + dr^, then the metric can be re-written 



d-\ 


-f{r)dY + ^ {dxif 
i=\ 


_dY_ 
Yf{r) ’ 


(3.4) 


which is in the same form of Eq. (4.1) of Ref. [33]. Sometimes 
it is more convenient to work with a compactified radial coor¬ 
dinate, and normalize the boundary coordinates by the scale 


of the black brane horizon. With u = I/{brY, t = t%%T/d, 
x‘ = x' %nT jd and /(«) = 1 — u‘^1^, the metric becomes 



{AnT/df 


d-\ 


—f{u)dY + ^ {dx^Y 


=1 


-b 


1 

4m 


d-\ 

—f{u)dY -b {d^Y 
i=l 


-b 


4m2/(m) 


AYf{h 


dY. 


dY 

(3.5) 


To derive the gravity/fluid correspondence, we take as our 
starting point the uniformly boosted black brane (3.1). 


B. Gravity/fluid correspondence 

To each asymptotically AdS bulk solution there is an asso¬ 
ciated metric and conserved stress-energy tensor on the time¬ 
like boundary of the spacetime at r —oo (see, e.g.. Ref. [34]). 
The boundary metric, in the case of (3.1) is Ij^v, while the 
boundary stress-energy tensor is 

This describes a perfect fluid with energy density p and pres¬ 
sure p given by 


d-l 

(3.7) 

^ l67tGd+ibd' 

1 

lenGd+ib^' 

(3.8) 


The stress-energy tensor is traceless, with equation of state 


P = 


P 

d-l' 


(3.9) 


as required by conformal invariance. Imposing the first law of 
thermodynamics, dp = T di, as well as the relation p+p = sT, 
gives the entropy density s and fluid temperature T, 


s=AT‘^-\ (3.10) 

P = ^AtY (3.11) 

d 

Here, A is a constant of integration. This is fixed to A = 
{AnY/ (167rGrf+i£/‘^^*) by equating T with the Hawking tem¬ 
perature. 

At this point, the fluid we have described is of constant den¬ 
sity, pressure and velocity. To go beyond the uniform fluid, b 
and are promoted to functions of the boundary coordinates 
xb. Importantly, these will be assumed to vary slowly; that is, 
if L is the typical length scale of variation of these fields, then 
b. With non-constant boundary fields, the metric (3.1) 
no longer describes a solution to the Einstein equation. How¬ 
ever, a solution can be obtained by systematically correcting 
the metric order by order though a derivative expansion, so 
that the Einstein equation is solved to any desired order in 
derivatives. One can then compute the boundary stress-energy 
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tensor corresponding to the metric at each order, and take this 
as defining the boundary fluid. 

After a rather long, but direct, calculation, the resulting 
boundary stress-energy tensor (to second order in derivatives) 
is 

^ (du^uy + ri^y) +n^v, (3.12) 

where the viscous part Il^v is (see, e.g., Eq. (3.11) of Ref. [4]) 
n^v = ~2T]CT^y 

+ 277Tn daOfiv) -f 

+ (AlC7^aCTv“-f A2C7;iaCUv“+A3a)^aCUv“)- (3.13) 

The shear and vorticity tensors are defined as, 

jiv = {d^ lUv) , (3.14) 

^^v = Pii Pv^ (3.15) 

We have employed angled brackets to denote the symmetric 
traceless part of the projection orthogonal to m^. 



(3.16) 

and defined to be the spatial projector orthogonal to u^, 

Ppv = Vpv + Ufj^Uy. 

(3.17) 

Notice that fl^y is symmetric and satisfies 



(3.18) 

= 0 . 

(3.19) 


The transport coefficients {tj, Tn, A,} for various dimensions 
can be found in, e.g., [6, 35, 36]. In particular, Vj = s/ (4;r). 

Projection of the Einstein equation along the boundary di¬ 
rections shows that the boundary stress-energy tensor is con¬ 
served, giving rise to the fluid equations of motion, 

0 = u^dyp + -^pdyu''-u^d'^Yluv, (3.20) 

d—\ 

(3.21) 

d—l 

The gravity/fluid correspondence thus provides an explicit 
link between black hole perturbations in the sufficiently long 
wavelength regime —described by small wave numbers—and 
relativistic hydrodynamics. Ordinary perturbation theory, by 
contrast, provides a solution that is valid for sufficiently small 
amplitudes, but cannot easily capture the transfer of energy 
between modes. Our coupled-oscillator approach in contrast 
does capture the leading mode-mode couplings that are mani¬ 
fest in the fluid picture, and it is in that sense valid for larger 
amplitudes (see Sec. IIB). As illustrated in Eig. 1, there is an 
overlapping regime where the predictions of both approaches 
can be compared. 



Figure 1. An illustration of the hydrodynamical expansion (small 
wave number) and black-hole perturbation (small amplitude). They 
both admit effective coupled oscillator descriptions. In AdS black- 
brane spacetime we compare the results from both sides of the dual¬ 
ity, in the shaded region of the plot. For small perturbation amplitude, 
this comparison has been done in the linearized perturbation theory 
(as depicted by region “A” and see for example [33]). For larger 
perturbation amplitude (region “B”), we are able to expand the com¬ 
parison to equations of motion with nonlinear couplings using the 
coupled-oscillator model. 


C. Mode expansion of the boundary fluid 

We now proceed to re-write the fluid equations as a set of 
coupled oscillator equations so that they can be compared with 
the equations we will derive on the gravity side. We denote 
the four velocity u^ = (7, u), where 7^ = 1 -Fu • u, and the 
density p = poe^. Keeping viscous terms to linear order in 
u and ^, and inviscid terms to quadratic order (as needed for 
the comparison), the energy conservation and Euler equations 
reduce to 

O^dr^+u-V^ + ^idrY+V-u), (3.22) 

d—l 

0 = + Vu+ 

--f'^v2u-F^V(V.u)) . (3.23) 

Po \ d d J 

Eurthermore, dropping nonlinear terms, 

0 = + (3.24) 

d—l 

d 

f^^vV')-F^^V(V.u('))V (3.25) 
Po \ d d J 

Linearized solutions are decomposed into two families of 
modes: sound and shear. A sound wave of momentum k takes 
the form 

^Ah{k)e-“^Pe‘'^-^k, ~ Bi,(k)e-‘^'’'e ''^-^. (3.26) 
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By solving the linearized equations (3.24) and (3.25), the dis¬ 
persion relation is found to be 

= + (3.27) 

Vfl — 1 “ Po 

and 

Bb{^) = ^—Ab{y^). (3.28) 

a —\ (Oh 

For the shear modes, ^ =0 and 

(3.29) 

with Ui • k = 0. The resulting dispersion relation is 

(0, = -i‘^^—k^ + 0{k^), (3.30) 

d po 

so shear modes are purely decaying. The general solution to 
the linearized fluid equations is simply a sum over sound and 
shear modes of different k and shear polarizations s. 

We are now in a position to include the effects of nonlinear 
coupling terms. To do so, we express ^ and u as sums over 
linear modes, but we allow for the coefficients A and B to be 
functions of time. The velocity ansatz then takes the form 


fact, if we ignore all the sound modes in the relativistic hydro 
equation, the resulting Eq. (3.32) is the same as the one for 
incompressible fluid (Appendix B), and they share the same 
conservation laws in the Fourier domain. 

Equation (3.32) expresses the fluid as a collection of cou¬ 
pled oscillators, to be compared with (2.28) on the gravity 
side. In the next section we shall apply the general formal¬ 
ism of Sec. II to the AdS black brane spacetime and directly 
match its mode coupling coefficients (for the fundamental hy¬ 
dro shear quasinormal modes) to the shear-shear mode cou¬ 
pling coefficients in Eq. (3.32). One can apply the same proce¬ 
dure to verify the correspondence in the sound channel (which 
we have not written down). We will only address the shear 
modes, as the main purpose of this work is to formulate the 
coupled oscillator model and to illustrate its technical details, 
rather than to provide a full verification of the gravity/fluid 
correspondence. We envisage that this framework shall prove 
its unique value when studying gravitational interactions in 
spacetime without a clear gravity/fluid correspondence, or in 
cases where the hydrodynamical (long-wavelength) approxi¬ 
mation becomes too restrictive. 


IV. LINEAR AND NONLINEAR GRAVITATIONAL 
PERTURBATIONS OF THE AdSj BLACK-BRANE 


u(x,f) 


L 


q'6(k,f)k + ^q'i(k,f)u^ 

S 


e 


ikx 


(3.31) 


where q'j(k,f) = As(k,f)e^‘®'^ and qb{Kt) = A/,(k,f)e^'“'’'. 
The coefficients are of course subject to a reality condition. 
Inserting this expansion into Eq. (3.23), and projecting it onto 
a particular shear mode, we obtain 

5A(k,f) (3.32) 

= i E [w'(pT)-q][w(k,f)-W"(q,f)Rv'(pTMy'(q,f) 

p+q=k,R,A" 

+ E E {■■■)qb{i?,t)qb{q,t), 

p+q=k, s' p+q=k 

Notice that the left hand side has been reduced to simply the 
time derivative of A^ because the mode function satisfies the 
linearized equation of motion. The right hand side describes 
the nonlinear coupling between modes. 

The second and third terms (coupling coefficients unspec¬ 
ified) in Eq. (3.32) describe the mixing between the sound 
modes and the shear modes, as well as between two sound 
modes. The coefficients to these terms contain fast [exp(/a)f) 
type] oscillatory time-dependent factors, so their effects tend 
to average to zero during the longer time scales in which we 
examine the growth and decay of modes. On the other hand, 
the first term describes the mixing between two shear modes, 
and it trivially satisfies the “resonant condition” in the time- 
domain since 91 ( 0 ^) = 0. This results in significant energy 
transfer between shear modes (and had we been performing an 
ordinary perturbative expansion would have resulted in sec¬ 
ular growth). It is then natural to expect that the effect of 
sound modes is sub-dominant in the turbulent process of con¬ 
formal fluids, where the viscous damping is less important. In 


In this section we study gravitational perturbations about an 
asymptotically AdS black brane within the context of the cou¬ 
pled oscillator model. We adopt this particular example for 
two reasons: On the one hand, the boundary metric of the 
background spacetime is flat, which simplifies calculations 
when performing wave function projections. On the other 
hand, the gravity/fluid correspondence is well established in 
this spacetime, and this allows us to compare results obtained 
in the gravity and dual fluid pictures, as depicted in Eig. 1. 
In particular, we shall focus on the analysis of shear modes 
at both linear and nonlinear levels. We also fix the spacetime 
dimension to t/-|-1 = 5, although it is straightforward to gener¬ 
alize the analysis below to other dimensions. Eor calculations 
within this section, we make further simplifications by scal¬ 
ing the coordinates such that b = I, so the horizon is located 
at r = 1. This means that we effectively choose T = \/n so 
[see above Eq. (3.5)] 

= ki = 2ki. (4.1) 


A. Linear perturbation 

Linear quasinormal mode perturbations of AdS black 
branes have been thoroughly analyzed in [33]. There, the fun¬ 
damental (slowly decaying) quasinormal modes of the space- 
time were shown to be the same as the hydrodynamical modes 
of the boundary fluid. The analysis was performed using the 
coordinate system of Eq. (3.5), whereas for our purposes it is 
more convenient to use the ingoing coordinates of Eq. (3.3). 
As discussed in Appendix C, choosing different coordinates 
leads to different definitions for the modes. At the linear level 
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there exists a clean one-to-one mapping of modes in different 
bases as each quasinormal mode is a solution to the linear Ein¬ 
stein equation. However, when studying nonlinear perturba¬ 
tions, their projection with respect to a mode-basis associated 
to a different coordinate system leads to an expansion with 
a less direct identification. In Appendix C we illustrate this 
point with a simple example describing a scalar field propa¬ 
gating on Minkowski spacetime. 

As demonstrated in [33], linear perturbations of the AdS 
black brane can be classified into shear, sound and scalar sec¬ 
tors. In addition, as the boundary metric is flat, it is straight¬ 
forward to Fourier transform the metric components along the 
boundary coordinates. The same logic applies when we adopt 
ingoing coordinates. Without loss of generality, we consider 
a mode whose boundary-coordinate dependence is For 
shear perturbations, the relevant metric components are then 
hra,hvaTh^oi, where a —x,y. Without loss of generality, we 
choose the polarization a = x, and impose the radial gauge 
condition = 0, with M = {r,v^z,x,y}. Defining the auxil¬ 
iary variables 


-ife 


Hvc — ^zx 9 


Hvr = hv 


(4.2) 


the independent components of the linearized Einstein equa¬ 
tion take the form 


0 = 5r- 


dH, 


■ +ik 


dH^ 

+ r 


0 = k^H„-5r^f 


dr 

dH, 


dr^ 


(4.3) 


^ dr^ dv ' dvdr 


We can further simplify this system by defining the master 
variable, ‘F = SrH^x- This satisfies the master equation. 


- -p (SrV'P)' + -f 7r'F -P = 0, (4.4) 


where in this section we will often denote partial deriva¬ 
tives as (•)' = dr and (•) = d^. To look for quasinormal 
modes, we hrst take advantage of the time translation sym¬ 
metry of the equation to impose a time dependence (so 
'F — —ico'¥). Solving the remaining spatial equation with ap¬ 
propriate boundary conditions at the horizon and spatial in¬ 
finity gives rise to a set of quasinormal modes in the ingoing 
coordinates, and the frequency spectrum co{k). 

To analyze the horizon boundary, we multiply Eq. (4.4) by 
/ and take the horizon limit r —1. The wave equation be¬ 
comes 


{dl+2drdrJ^> = 0, (4.5) 

with two independent solutions, 

5,,'F = 0, and (5^,-p2(9v)‘F = 0. (4.6) 

The ingoing boundary condition for the quasinormal modes 
selects 


As r —oo we impose a reflecting boundary condition (since 
the spacetime is asymptotically AdS), so the metric perturba¬ 
tion is required to vanish. This means that we should at least 
expect h — 0{l/r) and 'F = 0{l/r^). 

The above discussion applies to all quasinormal modes of 
our system. However, the dual fluid captures only the longest 
lived shear and sound modes, which have O) —> 0 as k —> 0 
(known as the “hydro” modes). In order to compare our re¬ 
sults with the fluid we therefore restrict to k <C 1. We can then 
construct the eigenfunctions perturbatively in k (and (o). In 
this expansion, the leading order part of equation (4.4) is 

{Sr^fY'i' + 5/r^'F' -P + (r'^/j'^F' = 0. (4.8) 


After imposing the horizon boundary condition, the solution 


'Fo = 


C(v) 


(4.9) 


where the subscript 0 indicates that this solves the leading or¬ 
der equation. (Notice that this solution also falls off suffi¬ 
ciently rapidly at spatial infinity.) To look for quasinormal 
mode solutions we take C(v) = 

The leading order solution 'Fq then sources the first order 
correction ‘F i through 

(5rV)'‘Pi + 5fr^% + rV'P'/ + 

= -{-k^^O + 7r'¥o + 2r^%). (4.10) 

The combined solution 'F = ‘Fq -P 'Fi is then 

1 -(k2 - 4/0)) log( 1 - r) - (k2 -p4/0)) log( 1 -p r) 


‘F = 




16r^ 

8 / 0 )arctanr-Pk^log(l -Pr^) 


16r^ 


(4.11) 


In order to satisfy the horizon boundary condition we must 
impose k^ = 4/0), resulting in 

2k^arctanr —2k^log(l -Pr) -Pk^log(l -Pr^)’ 

fj 


^> = 


^5 

(4.12) 

Using Eq. (4.1), we verify that k^ = 4/0) is equivalent to 
(a = —/k^/ 2 , which is exactly the dispersion relation of shear 
hydro quasinormal modes derived in [33] using a different co¬ 
ordinate system. In addition, it is easy to check that the dis¬ 
persion relation matches (3.30), derived on the fluid side. 

Knowing 'F, it is straightforward to use Eq. (4.3) to recon¬ 
struct the metric perturbations. For the shear modes consid¬ 
ered here, the metric perturbation is 


/, -_ ^ -imv+ikz 

““ 4^2 


1 + ——{nr -4r + 2) 
16 


16 


-(r^ — 1) ( 2 arctan r -P log 


l-Pr^ 

(TT)f 


(4.13) 


hvi = i—kre 


2„—i(Ov+ikz 


1 arctan r 


n 

4~Y~~ 

2\n I 


1 (i+PKl+p 

4 n 


-)■ 0 , r —)• 1 . 

dr^ 


(4.7) 


2 


(4.14) 
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B. Mode projection 


Having carried out the linear analysis, we are almost ready 
to calculate the shear-shear mode coupling coefficient. There 
is one more problem to tackle however, which is to project 
the Einstein equation onto an individual mode to see how a 
source term affects its evolution. As described in Sec. IIC, 
we adopt a technique that has been proven very powerful in 
solving similar problems [7, 25-28, 37]. Namely, we enlist a 
suitable bilinear form to project the equation onto individual 
modes. 

For later convenience, we define (j) = so that Eq. (4.4) 
takes the form 


Fourier transforming the wave operator in v, we define 






7 io) ^ 

—j-(i) —2i(0r 
n 



/ 


(4.16) 


We also define a generalized inner product, 

{x\n) = J^ drxn- 


(4.17) 


The operator Ha is not symmetric under this bilinear form, 
i.e., {x\Ha'if]) 7 ^ (//(BXlq), because of the fourth term in Ha- 
However, in the hydrodynamic limit (k ^ 1) this term is ne¬ 
glected, so (4.17) is suitable for our purpose of comparing to 
the dual fluid. 

For completeness, we note that should the need arises for 
the study of perturbations of higher overtones away from the 
hydro limit, we may use an alternative bilinear form (depen¬ 
dent on co) with respect to which Ha is symmetric so that 
{x\H(ori)a = {HaX\l)a>- In this case. 


= ^ drga{r)xri, with (4.18) 


2 _ f 

loggfi,(r) = -i(0 ( arctanr-f - log - 


+ r 


- const, 


is the unique option. There is one ga for each co, so we have a 
family of such generalized inner products. Using ga to project 
onto the mode with frequency CO [followed by a diagonaliza- 
tion procedure as per the discussion above Eq. (2.17)] is a nat¬ 
ural choice, and indeed leads to agreement with the Green’s 
function method for projecting modes (see Appendix D). In 
any case, to 0[k^), these generalized inner products reduce to 
Eq. 4.17. 

For the purpose of the time-domain analysis in the next sec¬ 
tion, we expect the effect of non-hydrodynamical modes [see 
Eq. (4.26) below] and the excitation of residual parts^*^ to be 


'* The prompt piece of the residual can be intuitively understood as the source 
terms propagating on the light-cone. Also notice that the source terms, as 
represented by Eq. (3.32) or Eq. (4.29), are linear in the hydrodynamical 
momentum, so overall the source terms are of 0(k), as is the excitation 
amount of the prompt residual. 


at least 0{k). Therefore only the hydrodynamical modes are 
important to leading order and we shall adopt the generalized 
inner product (4.17) for calculations, as it is easier to imple¬ 
ment in the time-domain analysis. As an example, we show 
below that this inner product generates the correct leading or¬ 
der (in k) frequency in the eigenvalue analysis. 

Let us now consider a simple example that demonstrates 
the essence of how to utilize this inner product to carry out 
perturbation studies. Suppose we perturb k to k-fe5k (e ^ 1) 
and ask for the change of co. On the one hand, based on the 
dispersion relation co = —ik^/4, we immediately know that 
5co = —ik5k/2. On the other hand, we can arrive at the same 
conclusion through a perturbation analysis of the eigenvalue 
problem defined by Eq. (4.16). 

The change k —7 k -f e5k causes H to pick up an extra 
term, —2ek5k/r^. We expect both the eigenfrequency and 
the eigenfunction to also change to order e, 

(p —>■ (j) +0{e ^), 

co^ co + e5co + 0{e^). (4.19) 


Plugging into the wave equation Eq. (4.16), and projecting 
both sides onto (j) while keeping only the (9(e) terms, we can 
eliminate the unknown function 0^*) to obtain 


—i5co —2k5k' 


\l/r^) 


{^\1 ^ + 2r^ / r^y 


-0(e) 




(4.20) 


which is consistent with our expectation. We note that it was 
necessary in this analysis to use the symmetry property of Ha 
to eliminate terms involving Although somewhat exces¬ 
sive for this simple problem, we see that with the help of our 
generalized inner product, it is now possible to carry out a per¬ 
turbation analysis in a manner analogous to the application of 
perturbation theory in quantum mechanics [38] (for a direct 
mapping of a wave equation with outgoing boundary condi¬ 
tion into a Schrodinger equation with non-Hermitian Hamil¬ 
tonian, see [24]). 


C. Nonlinear analysis 


We are now in a position to move beyond the linear 
level and study the second order (nonlinear) Einstein equa¬ 
tion (2.10). We begin by considering its projection onto the 
shear sector with spatial dependence and spatial polariza¬ 
tion a = z (see Sec. IV A). (It is straightforward to perform 
this projection onto a Fourier basis element with an ordinary 
inner product. The nontrivial aspect is the subsequent pro¬ 
jection onto the hydro mode.) The non-vanishing vx and rx 
components of the Einstein equation take the form 


5r- 


-hik^--hr = Tr; 


k^H,,-5Pf 


— Tv: 


dr 

dH, 


dr^ 


vz 4 ^VZ 


-[- ik - 


dH-, 


zx _ ^2 d^dlvz 


dvdr 

(4.21) 
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Figure 2. An illustration of three wave numbers satisfying the “mo¬ 
mentum matching” condition. 


We have formally written the nonlinear terms as “sources” on 
the right hand side of the equation. At quadratic order the 
nonlinear terms are 

Xrz = - T,, = -. (4.22) 


The inner product (•, •) is the ordinary inner product over the 
boundary spatial coordinates. Equation (4.21) is simply (4.3) 
with nonlinear terms included, and a simple switch of coordi¬ 
nates z. 

Since the second order Ricci tensor is a quadratic func¬ 
tion of the metric perturbation, which can be expanded over 
Fourier modes (and scalar, sound, shear sectors), the projec¬ 
tion (4.22) enforces a wave number matching condition on 
the terms that can contribute to the right hand side of (4.21). 
Namely, modes with wave numbers p and q can only act as a 
source for mode k if p-f q = k (see Fig. 2). [This of course 
also holds for the fluid analysis in (3.32).] We define the an¬ 
gles 01 = arccos(^ • k) and 02 = arccos(p • k). 

Following the same procedure as in the linear analysis, 
we re-write Eq. (4.21) in the form of a sourced version of 
Eq. (4.15), 








— frz 


(4.23) 


Since only first order time derivatives appear in this wave 
equation and we know from the previous subsection that the 
quasinormal frequency is purely imaginary, this shear hydro- 
dynamic mode belongs to the class described by Eq. (2.29). 

We now proceed to compute the nonlinear source Sj 
[see (2.29)] using the generalized inner product of Sec. IV B. 
As in Sec. 11C, we first express the field 0 as a sum over radial 
modes 


0 ^AQ{t)e-'^'’xo{r) + Y^Aj{t)e-'‘^A^.^r)+o{k) , (4.24) 
j>o 


accordingly The 0{k) appearing in the expression for 0 
includes the residual contribution under the hydrodynamical 
approximation (see Footnote 4). 

We can now plug (4.24) into the wave equation (4.23), and 
then take the generalized inner product of both sides with Xo 
using (4.17). Within this computation, the effect of the non- 
hydrodynamical terms is at least 0{k) [in fact 0{k^)] as we 
claimed in Sec. IV B. This is because solves the lin¬ 

ear equation (4.15), so for j >0 



= 0{k^). (4.26) 

Given this observation, it is now simple to show that 

4 {Xo\'^Xo/r^+ '^r^{Xo/r^)') 

--\^vAr=x, (4.27) 

where we dropped high order [0[k^) and higher] terms in k, 
including nonlinear terms containing time derivatives (as dis¬ 
cussed in Sec. II C). Using Eq. (2.4), the mode expansion of 
h^v, and after some lengthy but nevertheless straightforward 
calculations, one can show that the shear-shear mode coupling 
coefficient arising from (4.27) is 

Kkpci = iksm{e2-9i), (4.28) 

which agrees with the result obtained with its fluid counterpart 
from Eq. (3.32) 

Kkpy = f[u,/ (p, f) • q] [uj(k, f) • u,.// (q, f)] + (p o q). (4.29) 

We end this section by noting that the agreement between 
the mode coupling coefficients inferred from the fluid equa¬ 
tions and the AdS black brane perturbation theory relies on 
the fact that they are computed using the same mode basis, and 
that the comparison is made in the regime where k I and 
\!i\ <C 1 (cf. Fig. 1). However, the coupled oscillator model is 
applicable more broadly. 


V. CONCLUSIONS 


but we allow for the modes to have additional time depen¬ 
dence through the mode amplitudes Aj. Here the spatial 
wavefunctions are denoted Xj, with j — 0 corresponding to 
the hydro mode. The non-hydro modes all have frequencies 
COj = (9(1), while cOq = —ik^/4. While 0 is to be matched to 
modes of the 4-velocity on the fluid side, we normalize the 
wave function 

Zo = 4 (4.25) 


The study of nonlinear wave phenomena is undoubtedly 
a fascinating subject. Gaining understanding in the partic¬ 
ular case of general relativity poses unique challenges even 


^ This is of course just an inconsequential overall constant rescaling of A, the 
more important goal is to match the angular dependence of the coupling 
constants. 
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given the fixed speed of propagation of physical perturbations. 
These challenges are rooted in the covariant nature of the the¬ 
ory and physical degrees of freedom often hidden within a 
larger set of (metric) variables. These issues have hampered 
understanding of gravitational perturbations beyond linear or¬ 
der except in a few specialized regimes [12-16], seamingly 
leaving full numerical simulations as the main tool to try to 
understand these issues (for a recent overview of these efforts, 
see [39] and references cited therein). 

In the current work, we have presented a model to capture 
the nonlinear behavior of gravitational perturbations®. This 
model regards the system as composed of a collection of non- 
linearly coupled (damped) harmonic oscillators with charac¬ 
teristic (isolated) frequencies given by quasinotmal modes. 
By construction this model reproduces standard results ob¬ 
tained at the linearized level. At the nonlinear level, it de¬ 
scribes mode-mode couplings and their effect on frequency 
and amplitude shifts. As an illustration, we have shown 
how our model reproduces recent results captured through the 
gravity/fluid correspondence via a purely gravitational calcu¬ 
lation. Importantly, the applicability of our formalism is not 
restricted to long-wavelength perturbations—as in the case 
of the gravity/fluid correspondence—so the coupled oscilla¬ 
tor model can also treat so-called “fast (non-hydrodynamical) 
modes” of perturbed black holes [40]. As a consequence it 
can be employed to study a broder phenomenology than that 
reachable via the correspondence^. We stress that our for¬ 
malism is also applicable beyond asymptotically AdS space- 
times. Thus it can also help shed light on nonlinear mode 
generation in perturbations of asymptotically flat black hole 
spacetimes [42, 43]. 


Appendix A: Brief overview of coupled oscillator systems 

Consider a family of nonlinearly coupled harmonic oscilla¬ 
tors governed by. 


= Y^C^^jlIqkqi + >^^^iqkqi + >^fkiqkqi) = Sj , (Ai) 

kl 

where (ojqj is the restoring force and jj is the damping coef¬ 
ficient. Each oscillator’s displacement can be decomposed in 
the same way as Eq. (2.14), with coj satisfying 

- coj - ijjCOj + 6}j = 0. (A2) 

In the presence of nonlinear mode-mode coupling ^ 
0), Aj and Bj are both time-dependent. In fact, we can take 
one more time derivative of the first equation in Eq. (2.27) , 
and obtain 


{Aj - i(OjAj)e qj + iqj) 

UJj -t- UJj 

= + f 4- + 

C0j + 0)J 1 ^ 0 ); + ®]' i(o)J+COj) J ’ 


COj + CO* 


i(0j 

COj + CO* 


{(0*qj + iqj), 


(A3) 


such that 
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and similarly 


Bj = - 


iSj 


COi + CO 


-m-r 
re 1 . 


J 


(A5) 


These effective equations of motion have the same kind of 
first-order form as Eq. (2.29) and Eq. (2.28), which means 
that one can utilize results from previous studies on nonlinear 
coupled oscillators to analyze nonlinear gravitational interac¬ 
tions. 


® In this work we have included up to three-mode interactions, but the for¬ 
malism can be extended to include higher order interactions. 

^ Recently, resummation techniques have been proposed to take some of 
these higher modes into account within an extended hydrodynamical de¬ 
scription [41], This requires knowledge of the hydrodynamical expansion 
to very large orders. 


Appendix B: Two-dimensional incompressible flnid in the 
inertial regime 

Here we review the Navier-Stokes equation for a two- 
dimensional incompressible fluid. This discussion highlights 
how a new symmetry for the mode-mode coupling coefficient 
arises in the mode-expansion picture. Such symmetry is crit¬ 
ical for the double-cascading (inverse energy and direct en- 
strophy cascades) behavior in two-dimensionalfluids. A more 
detailed discussion can be found in Ref. [44]. 
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The Navier-Stokes equation for an incompressible fluid in 
the spatial-frequency domain reads 

(^ + Vk2) uj{k,t) 

= ikiPj„{k) Y, u„{p,t)ui{q,t) 

p+q=k 

ikjP (k) 

=—w - E {'^nip,t)ui{q,t) + u„{q,t)u[{p,t)] (Bl) 

p+q=k 


where u(x,f) = ^i 5 e'’^’‘u(k,r) and Pjn{k) = 5j„ — kjk„/k^. In 
incompressible fluids, the condition V • u = 0 translates to k • 
u(k,f) = 0 in the Fourier domain. We can write u(k,f) as 

u(k,f) = A(k,f)M(k,f) , (B2) 

where u{k,t) satisfies ii k = 0 and m-m = 1. In 2 +1 fluids, u is 
unique for any k. Using the new variables, the Navier-Stokes 
equation can be rewritten as 


^ + vk^]A[k,t)=i Y «■(■«■,P,q)A(p,f)A(q,f) 
/ p+q=k 


= i Y {{u{Kt)-uip,t)][^-uiq,t)] + [u{k,t)-u{q,t)][k-u{p,t)]}A{p,t)A{q,t). (B3) 

p+q=k 


This is the same as the shear-shear coupling term in Eq. (3.32), which is already written in a form consistent with the coupled 
oscillator model. 

In the inertial regime we shall set the viscosity coefficient v to zero (as such coefficient only governs the extent of the regime 
but not the behavior within it) and recall that u(x,f) must be real. One can then show that 


d[uj{k,t)u*{k,t)] 

dt 


Y Im{[u(k,f) • u(p,t)][k • u(q,t)] -f [u(k,r) • u(q,f)][k • u(p,t)]} 

p+q+k=0 

= Y im[K'(k,p,qM(p,f)^(q,fM(k,f)]- 

p+q+k=0 


(B4) 


Energy conservation requires that 


d[uj{k,t)u*{k,t)\ 

dt 


d[uj{p,t)u*{p,t)\ d[uj{q,t)u*{q,t)] 


dt 


• + ' 


dt 


= 0 , 


(B5) 


which is equivalent to demanding 


K'(k,p,q) + K-(q,k,p) + K-(p,q,k) =0 (B6) 

for any vectors k, p, and q satisfying p + q -f k = 0. It is straightforward to check that the above relation is automatically 
satisfied given the expression of Moreover, for 2 -f 1 fluids, by using the fact that k • u(k, f ) = p • u(p, f ) = q • u(q, f) =0 and 
the identity 

sin^ 01 cos(02 — 03 ) + sin^ 02 cos (03 — 0i) -l-sin^ 03 cos( 0 i — 82 ) = 0, (B7) 


for V 01 -f 02 + 03 = we can show that an additional sym¬ 
metry for the mode-mode coupling exists, which is 

k2K-(k,p,q)-fg'2K-(q,k,p)-f/K-(p,q,k) =0. (B8) 

This additional symmetry is directly connected with the ad¬ 
ditional conserved quantity in 2 -f 1 fluids: enstrophy. With 
two conserved quantities in the inertial regime, Kraichnan 
[44] explained that a dual-cascading behavior should be ex¬ 
pected in the turbulent regime. This example strongly sug¬ 
gests that the symmetry of the mode-mode coupling coeffi¬ 
cients in our coupled oscillator model could be crucial for 
classifying the nonlinear behavior of gravitational evolutions. 


Appendix C: Expansion in two different bases 

Let us imagine a simple example of a scalar held whose per¬ 
turbations propagate on a 2-dimensional flat spacetime with 
time-like boundaries at x = 0 and x = 1. Eor comparison pur¬ 
poses, we have assigned two coordinate systems in this space- 
time: standard Cartesian coordinates (f,x) and “null” (v,x) 
coordinates, with v = f -|-x. Eor simplicity, we impose Dirich- 
let boundary conditions 't>|x=o = = 0 for the wave. At 

linear order, the scalar wave satisfies the following wave equa¬ 
tion 

(Cl) 


(_52 + 52 )^ = 0 
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'' constant v 

\ 


constant t 


x=0 x=l 

Figure 3. An illustration for mode decompositions of a scalar field in 
a flat spacetime. At each point (such as the star in the diagram), we 
show two possible mode bases with respect to which to decompose 
the scalar wave. 


in the (v,x) coordinates. We note that the mode amplitudes 
are generically time-dependent now. 

Pick an arbitrary point in the spacetime (for example, the 
one labeled with a “star” in Fig. 3). There we can ask whether 
the matching described in Eq. (C5) still holds for the two dif¬ 
ferent mode expansions at that point. As we can see from 
Fig 3, these two mode expansions sample two different slices 
of the spacetime: one at constant f and the other at constant v. 
Unlike the linear case, the scalar wave distributions on these 
two slices can be made quite “independent” of each other by 
freely detuning the nonlinear terms in the wave equations. In 
the end, the largely independent data on these two slices im¬ 
ply that simple mappings such as Eq. (C5) no longer exist for 
mode expansions under different bases in the general nonlin¬ 
ear scenario. However, we emphasize that despite the lack 
of a simple mapping between them, both mode expansions 
are equally valid in describing the wave evolution. Although 
our present analysis is performed using this simple example 
where the mode expansion is complete, we see no reason why 
a similar conclusion would not hold for quasinormal mode ex¬ 
pansions of generic spacetimes. 


in the {t,x) coordinate system or 

{d^ + 2d,d,)<P = 0, (C2) 

in the (v,x) coordinate system. 

Based on the wave equation and the boundary conditions, 
we can see that this is a standard Sturm-Liouville problem, 
where it is straightforward to write down the solutions of the 
wave equation in a mode expansion 

^{t,x) = ^ +Bje‘^J‘) sm{jKx) , (C3) 

J 

and 

d>(v,;r) = ^ sinijnx ), (C4) 

J 

with cOj = jn. It is obvious that we can match up the linear 
modes from the two different expansions above, and in fact 
we can make the identifications 

Aj=Aj, Bj=Bj. (C5) 

Now suppose nonlinear terms (O^, <I>^ or even higher or¬ 
der) are present in the wave equations, resulting in a new so¬ 
lution [<I>(f,.«) or <I>(v,x)]. For such a wave, we can still choose 
constant-f or constant-v slices, and use the above spatial mode 
basis to perform a decomposition 

= E [Ajit)e^'“A sin(;7rx), (C6) 

J 

in the (f,x) coordinates, and 

4>(v,x) = E [A;(y)e-'“A-e'“A+B sin(;;rx) 

(C7) 


Appendix D: Coupled oscillator model in Schwarzschild 
spacetime 

As discussed in Sec. II, generic linear metric perturbations 
can be decomposed into quasinormal modes plus a resid¬ 
ual part. Unless we are dealing with normal modes which 
form a complete basis, or under certain physical conditions in 
which quasinormal modes dominate (e.g., AdS perturbations 
in the hydrodynamical limit), ignoring the contribution from 
the residual part should always require justification. Here we 
offer an alternative way of arriving at the coupled oscillator 
model, using the Green’s function approach (see also [9]). Us¬ 
ing this method, the quasinormal mode excitations can be un¬ 
ambiguously determined given a driving source term. So far 
this approach can only be demonstrated for perturbations with 
separable wave equations, such as Schwarzschild and Kerr 
perturbations, and we shall leave extensions to more general 
spacetimes to future studies. 

To simplify the problem, we assume that the angular de¬ 
pendence has been factored out, and we focus on the nonlin¬ 
ear evolution of modes with spherical harmonic indices {I, in), 
which satisfy the Regge-Wheeler (odd partity) and Zerilli- 
Moncrief (even parity) wave equations 

Here r* = r-|-2Mlog[r/(2M) — 1] and 'Fe,‘Po are the Zerelli- 
Moncrief and Regge-Wheeler gauge invariant quantities, re¬ 
spectively. The expressions for the potential and angular- 
projected source can be found in [45, 46]. In our present 
study, 5e/o is defined by the second order Ricci tensor, which 
is bilinear in the metric perturbations. 

Without the source term, for fixed time dependence 
there are two independent solutions to each wave equation. 


'I'e/o = 5e/o(6f)- (Dl) 
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One solution asymptotes to 


with 


U\n —y 6 


-(CD(r+r*) 


r* —>■ —oo 


(D2) 


near the event horizon, and 


Dn — 2CU|| 


dC{Y\ 

d(0 


^out(®n) • 

a>„ 


(D8) 


at spatial infinity. The other solution satisfies 


Mout ^ r*—>-00 (D4) 

at the spatial infinity, and 

Mout ^ CUco)e-‘‘^^‘+'-*^ +Couiico)e-'“^'-^^\ r, ^ — 

(D5) 

near the horizon. At the quasinormal mode frequencies 
(On, these two solutions become degenerate, and Cin(ft)n) = 

^OUt(^H) ~ 0- 

Using the Green’s function technique. Leaver [19] showed 
that ‘P can be decomposed as 


‘P = ‘Pqnm + ‘Pf + ‘Pbc, (D6 ) 


Notice that we are taking the real part because this QNM con¬ 
tribution is supposed to sum over both positive and negative 
frequencies. Also note that in order to maintain causality, we 
have introduced an upper bound t into the time integral of 
Eq. (D7), while in the original paper [19] this bound was set 
to oo (see also [47]). From Eq. (D7), it is then straightforward 
to derive the equations of motion for the amplitude of mode n 


^iCOn^ r 

An{r,t) = —— / dr'^Uin{r')S{r',t) = ——{uin\S)a., (D9) 
Un J 


where ‘Pp is the contribution from high-frequency propagator, 
‘Pbc is the branch-cut contribution in the Green function cal¬ 
culation, and 'Pqnm is the quasinormal mode contribution that 
we seek. In addition, he showed that 


‘Pqnm(tO = 2 Re 


E 


Min(r)e 


Dn 





(D7) 


where the integration should be performed as a contour inte¬ 
gral in the complex r' plane to ensure convergence [27]. In¬ 
terestingly, when we apply this Green’s function technique to 
analyze generation of the shear quasinormal modes in Sec. IV 
(as the wave equation is separable), we find that the gener¬ 
alized inner product (•|•)cI coincides with ('I')® defined in 
Eq. (4.18). 
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